Effects of nucleus initialization on event-by-event observables. 
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Abstract 

In this work we present a study of the influence of nucleus 
initializations on the event-by-event elliptic flow coefficient, 
V2. In most Monte-Carlo models, the initial positions of the 
nucleons in a nucleus are completely uncorrelated, which can 
lead to very high density regions. In a simple, yet more re- 
alistic model where overlapping of the nucleons is avoided, 
fluctuations in the initial conditions are reduced. However, V2 
distributions are not very sensitive to the initialization choice. 

1 Introduction 

The main goal of studying ultra-relativistic heavy ion re- 
actions is to identify (or not) the formation of the Quark- 
Gluon Plasma (QGP). Event-by-event (EbyE) analysis 
of global observables can significantly contribute in this 
direction [T]. One way to measure EbyE fluctuations 
is the statistical approach: each observable in an event 
should contain fluctuations and the distribution of such 
an observable can be characterized by its mean value and 
higher moments [2]. 

The most interesting fluctuations are dynamical, since 
they give important information about the formation of 
the system. Possible examples are the occurrence of jets, 
giving rise to fluctuations in the high pt tail of trans- 
verse momentum distributions, and the fluctuation in the 
anisotropic flow coefficient V2 due to unusual hard/soft 
equation of state or fluctuating initial conditions [TJ [3] . 

Other sources of fluctuations of statistical or technical 
nature can also be present in EbyE distributions. For 
example, finite multiplicity affects the determination of 
(pt), ratios of multiplicities of particle species and also 
the strength of the anisotropic flow coefhcient, due to 



' bemardo Sif. ufrj. br 



imprecise determination of event plane [T]. In this work, 
we are interested in a technical source of fluctuations, 
commonly present in event generators, the initialization 
of the positions of the nucleons in a nucleus before the 
collision. 



2 Hydrodynamical model 

We have proposed that dynamical and technical sources 
of fluctuations can be studied in the framework of rela- 
tivistic hydrodynamics, which allows one to separate the 
relevant physics from statistical noise [i]. If the local 
thermal equilibrium is attained in relativistic heavy ion 
collisions, hydrodynamical description may be most ade- 
quate for the space-time evolution of the system. Once 
the initial conditions (spatial configuration of 4-velocity 
field and conserved currents) of the system are specified, 
the principal factor which characterizes the hydrodynam- 
ical motion is the equation of state of the matter. When 
the hydrodynamical prescription is not valid anymore, 
one must employ a decoupling criterion to generate par- 
ticles, i.e, hadrons. So we have only three inputs in the 
hydrodynamical approach: initial conditions, equation of 
state and the decoupling criterion. 

In this work, we employ the numerical code SPHERIO 
[5] , which solves the hydrodynamical equations of motion 
in 3 dimensions and can deal with any kind of spatial 
configuration in the initial conditions. For the equation 
of state, we have adopted a quark-gluon free gas with a 
bag constant of 380 MeV/fm'^ for the QGP phase, and 
a hadron resonance gas with excluded volume (only for 
baryons) in the confined phase. Almost all the resonances 
up to 2.5 GeV have been used. The phase boundary is 
obtained via the Gibbs criterion, and the transition is first 
order for every chemical potential. As a decoupling crite- 



1 



rion, we adopted Cooper-Frye procedure [S] with a fixed 
freeze-out temperature of 140 MeV. For further details in 
the equation of state and decoupling prescription used, 
see [7]. 

3 Nucleus initialization: corre- 
lated versus uncorrelated nu- 
cleon positions. 

For the initial condition, wc use the NEXUS event gen- 
erator [5] . Based on the Gribov-Regge model of hadronic 
colhsions, it generates a spatial distribution of the energy- 
momentum tensor T^'^ and the baryon number density n b 
on the hyper-surface r = const. Monte-Carlo generation 
of events from this model give rise to physical fluctuations 
in the initial conditions. However, the procedure chosen 
to initialize the nuclei before the collision, may generate 
also unphysical fluctuations, which can influence EbyE 
observables. 

Usually one employs the Wood-Saxon distribution 



Pir) = Pa 



1 



l + exp[{r ~ R)/D]' 



(1) 



to determine the initial position of the A nucleons within 
a nucleus. The constants in cq. ([1]) are the nuclear density 
po = 0.16 fm~'^, the nucleus radius R and the diffuseness 
parameter D ^ 0.55 fm (usually). 

The nucleons determined this way are completely un- 
correlated, like a gas of free particles. So, it can hap- 
pen that several nucleons occupy the same position in 
space, which is unphysical. However, this method is used 
in many Monte-Carlo codes, since it is widely believed 
that unphysical fluctuations average out, at least for in- 
clusive spectra. Some models, e.g HIJING [8] are aware of 
this problem, and employ a rejection technique for nucle- 
ons which are closer to each other than a given distance. 
However, for realistic distances (in the order of 2 times 
the proton radius ^ 1.6 fm) this leads to a pushing of 
nucleon positions towards the shell and lowers density in 
the center. Nexus generates nuclei, with completely un- 
corelated positions. 

The most suitable initialization of a nucleus would be 
a selection of nucleon positions according to a realis- 
tic wave function which includes correlations [10| . Each 
nucleon would be surrounded by a corelation hole [TU] . 
which prevents other nucleons from occupying the same 
space. For simplicity, we employ a lattice model in or- 
der to demonstrate the effect. The nucleons are placed 
on a body-centered cubic lattice (BCC) and the posi- 
tions are then accepted with a probability p/po given 
by eq. H]). The packing efficiency of a BCC lattice is 



e = 8/3 * 7r/(4/V3)^ = 0.68 and the lattice spacmg IS 
I = (e//9o/7r*3/4)i/H4/V^= 2.32 fm [U]. This method 
correctly reproduces the density in the center of the nu- 
cleus, as can be seen in Fig. [T] To avoid lattice artifacts, 
each nucleus is rotated by some random angles. We shall 
call this method lattice (or correlated) initialization. 
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Figure 1 : Distribution of the nucleon density in the center 
of the nucleus (r < 4 fm) , for correlated (dashed line) and 
uncorrelated (solid line) nucleons. 



In Fig. [2]we show a plot for the energy density distribu- 
tion for the correlated (dashed line) and the uncorrelated 
(full line) cases. Each distribution represents 100 NEXUS 
events for Gold-Gold collisions with -y/s = 200 GeV per 
nucleon pair, with zero impact parameter b (central col- 
lisions). One can easily check that the energy density 
from correlated nucleons is narrower than the uncorre- 
lated case, as expected. It can also be seen that the aver- 
age energy density is a larger for the lattice type initial- 
ization. 

Au+Au@200 (b = fm, 100 events) 




energy density 



Figure 2: Energy density distribution for central collisions 
(b = 0). 
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4 Results 

In Figs. [21 m and we show EbyE distributions for the 
eUiptic flow coeflicient V2 (integrated over pt) of pions, 
computed for = 2 fm (central coflisions), 6 = 7 fm 
(mid-central) and 6 = 10 fm (peripheral collisions), at 
mid-rapidity, for Gold-Gold reactions at y/s ~ 200 GeV 
per nucleon pair. We show the histograms and the Gaus- 
sian fits for the observable originated from the correlated 
(dashed) and uncorrelated (solid) nucleus initialization. 
We computed 100 events per type of initialization at each 
impact parameter. Despite the differences in the energy 
density distributions, the V2 in all centralities shows very 
similar variances a. That may be caused by the hydrody- 
namical expansion which reduces the effects of the fluc- 
tuations. 

b =2fm, 100 NEXUS-ev, 1000 FO-ev (T = 140 MeV) 
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ELLIPTIC FLOW (all pions, y =0, pj < 3 GeV/c) 

Figure 3: Elliptic flow distributions for central collisions 
(& = 2 fm). 



b =1 fm, 1 00 NEXUS-ev, 1 000 FO-ev (T = 1 40 MeV) 




ELLIPTIC FLOW (all pions, y =0, pj < 3 GeV/c) 

Figure 5: Elliptic flow distributions for peripheral colli- 
sions (b = 10 fm). 

5 Conclusions and perspectives 

In this work we studied the effects of nucleus initializa- 
tion on EbyE observables. The V2 parameter is not very 
sensitive to the choice of nucleus initialization. The influ- 
ence on other observables are in progress. We also plan 
to study the effects of impact parameter distribution, re- 
action plane determination (for the V2 case) and other 
flnitc multiplicity related observables. The goal is to de- 
velop a method that separates physical from statistical 
fluctuations. 
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b =7 fm, 100 NEXUS-ev, 1000 FO-ev (T = 140 MeV) 
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ELLIPTIC FLOW (all pions, y =0, pj < 3 GeV/c) 

Figure 4: Elliptic flow distributions for mid-central colli- 
sions {b = 7 fm). 
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